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Abstract 

Quantum entanglement between qudits - the d-dimensional version of qubits - is rel- 
evant for advanced quantum information processing and provides deeper insights in the 
nature of quantum correlations. Encoding qudits in the frequency modes of photon pairs 
produced by continuous parametric down-conversion enables access to high-dimensional 
states. By shaping the energy spectrum of entangled photons, we demonstrate the creation, 
characterization and manipulation of entangled qudits with dimension up to 4. Their re- 
spective density matrices are reconstructed by quantum state tomography. For qubits and 
qutrits we additionally measured the dependency of a d-dimensional Bell parameter for var- 
ious degrees of entanglement. Our experiment demonstrates the ability to investigate the 
physics of high-dimensional frequency entangled qudit states which are of great importance 
for quantum information science. 

1 Introduction 

Entanglement is one of the most intriguing features of quantum theory It has been exper- 
imentally revealed by the observation of correlations with no classical origins. Specifically, 
through Bell inequalities, the non-locality of nature has been tested by numerous experiments 
using entangled two-dimensional states (qubits) [2J. Entangling d-dimensional states, denoted 
as qufiits, allows to formulate generalized Bell inequalities, which are more resistant to noise 
than their two-dimensional predecessors and lowers the threshold of the detection efficiency 
for loophole free Bell experiments [|3l SI |5l |6l. Entanglement is also a fundamental resource 
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for quantum information. Here, entanglement in higher dimensions is a prerequisite for imple- 
menting more complex protocols. For instance, the effective bit rate of quantum key distribution 
(QKD) can be enhanced, while at the same time the secret key rate and the robustness to errors 
increases with d [|3. 

Due to their low decoherence rate, photons are used in many experiments as a robust carrier of 
entanglement. Photonic entangled states are usually produced by the nonlinear interaction of 
spontaneous parametric down-conversion (SPDC) dSl. The coherence of this process, together 
with conservation rules, can generate entanglement in the finite Hilbert space of polarization 
states flU. Entanglement in infinite Hilbert spaces can be realized for transverse (momentum) 
or orbital angular momentum (0AM) modes [ITOl [TTl IT2l [TSl IT4l [T5l and for energy-time states 
[fT6l . The amount of entanglement is commonly quantified by the Schmidt number K. For 
transverse wave vector entanglement K is on the order of 10 for perfect SPDC phase matching 
conditions [fTTl . approximately 400 for specific non-perfect phase matching conditions ifTOl . 
and approximately 50 for 0AM entanglement [fTSl . 

Similar Schmidt numbers can be achieved in continuous energy-time entanglement generated by 
short pump pulses [|T6l[T8l but much larger K numbers are obtained for a quasi-monochromatic 
pump laser. In order to access these high-dimensional states, an experimental setup with many 
control parameters is required. In practice, the infinite Hilbert space is projected onto a finite 
space of, for example, discrete time- or frequency-bins. In the time-bin subspace with = 3, 4 
two-photon interferences have been observed by interferometers with multiple arms [[T9l l20l . 
This, however, requires interferometric stability and becomes prohibitively complex for higher 
dimensions. In the frequency-bin subspace, interferences between two entangled photons, each 
in an effective two-dimensional space, have been observed by manipulating the spectra with a 
combination of narrowband filters and electro-optic modulators [[2T|. In this approach a com- 
plex modulation scheme would be needed to address qurfits in higher dimensions. Besides 
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being very useful for QKD implementations, the aforementioned experiments do not provide 
sufficient control of the phase and amplitude of qurfits in the frequency domain to extensively 
study the properties of d-dimensional states with d > 2. 

Here, we demonstrate a methodology that allows for full control over entangled qurfits through 
coherent modulation of the photon spectra. It is derived from a classical pulse shaping ar- 
rangement and contains a spatial light modulator (SLM) as a reconfigurable modulation tool. 
This technique is widely used in ultrafast optics [|22] and has been adapted to manipulate the 
wavefunction of energy-time entangled two-photon states [|23l l24| . The flexibility of the ex- 
perimental setup enables the generation, characterization, and manipulation of rf-dimensional 
qudit states. We realize qudits with dimensions of up to 4. Specifically, we verify the genera- 
tion of the maximally entangled qurfit states through tomographic quantum state reconstruction. 
Subsequent Bell-type measurements demonstrate the applicability of the reconstructed quantum 
states. The versatility of the SLM based setup allows to test theoretical predictions beyond two- 
dimensional entangled states. As a first demonstration, we investigate the sensitivity of the Bell 
parameter for maximally and non-maximally entangled qubit and qutrit states by varying the 
degree of entanglement. 

2 Discretization of frequency-entangled photons 

We consider entangled photons generated in a SPDC process of type-0 where all involved pho- 
tons, the pump, the created idler (i) and signal (s) photon are identically polarized [[25l . Re- 
stricting the configuration of the three photons to the case where they are mutually coUinear, 
the entangled two-photon state can then be written as 




s 



(1) 
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where A{ui, Ug) = a{ui, Us)^{cOi,cOs) describes the joint spectral amplitude of SPDC in terms 
of the pump envelope function a{uji, ojs) and the phase matching function cOs). Idler and 
signal photons with corresponding relative frequency Ui^s are created by the operators a| ^(wj s), 
acting on the combined vacuum state |0)j|0)s. We have calculated the entropy of entanglement 
E = — Tr(plog2p) of A{uji,ujs) to be _E = (22.0 ± 0.2) ebits for a pump bandwidth of 5 
MHz by means of a numerical approximation method [|26l l27l . Here, p denotes the reduced 
density operator of the idler or signal photon subsystem. This amount of entropy is the same 
as in a maximally entangled quc/it state with d = 2^ ^ 4.2 ■ 10^. As a further quantification 
of entanglement, the Schmidt number K = 1 /Tr(/3^) has been computed numerically to be 
K ^ 2.4 ■ 10^. In order to use this large resource of entanglement for quantum information pro- 
cessing, we encode qurfits in the frequency domain by projecting the state \iJj) into a discrete (P- 
dimensional subspace spanned by the states \j)i\k)s with |j) . ^ = dufj"^ (w) al^^{u) |0) • ^ 
and j = 0, d — 1. The projected state is then expressed by 

d-l d-1 
j=0 k=0 

with coefficients Cjk = duidugfj* (uii) f^* {(^s) A (wj, Ug). The functions /j''^ (u) can 

be chosen arbitrarily under the condition to be orthogonal dcu fj'^* (u) fl''^ (cu) = 5jk- For 
the experiments presented here, we specifically define frequency-bins according to 



Imposing \ujj — uJk\ > {Auj + Auk) / 2 for all j, k ensures that adjacent bins do not overlap. For 
simplicity, we assume in the following a continuous wave pump by a{ui, Ug) oc 6{uji + Ug) and 
therefore \ip)^'^^ becomes restricted to its diagonal form 




1 / ^/Aujj for — ujj\ < Auj 
otherwise. 



(3) 



d-l 




(4) 



i=o 
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3 Experimental setup 



To experimentally generate type-0 entangled photons degenerated at 1064 nm, we pump a 11.5 
mm long positive uniaxial and periodically poled nonlinear KTiOP04 (PPKTP) crystal with a 
poling periodicity of 9 /im by means of a quasi-monochromatic Nd: YVO4 (Verdi) laser centered 
at 532 nm featuring a narrow spectral bandwidth of about 5 MHz (Fig.[il). The coUinear pump 
beam is focused into the middle of the PPKTP crystal with a power of 5 W. Mounted in a 
temperature stabilized copper block at 30.7 °C, the down-conversion crystal creates entangled 
photons with a spectral bandwidth of 105 nm. The measured photon flux power is 0.8 /iW and 
corresponds to a spectral mode density of 0. 15. This ensures that we are below the single photon 
limit and, therefore, the entangled photon pairs are temporally separated from each other [|28l . 
To compensate for group-velocity dispersion in the setup and to allow for coherent shaping of 
their spectra, idler and signal photon are imaged through a four-prism compressor arrangement, 
where the first prism deflects the residue of the pump into a beam dump. At the symmetry 
axis between the second and the third prism a SLM (Jenoptik, SLM-S640d) is aligned along 
the spatially dispersed down-converted spectrum. This device consists of two similar nematic 
liquid crystal arrays of 640 pixels, each with a width of 100 /im and separated by a gap of 3 
fim. The orientation of the liquid crystal molecules within a pixel can be controlled by the 
applied voltage. Together with a linearly polarized input beam and a polarization dependent 
detection scheme, the phase and amplitude of the transmitted frequencies at each pixel can be 
modulated [22J. The effect of the SLM on each photon is described by a complex transfer 
function M* ■^(w). A frequency-bin structure according to Eq. (|3) is then implemented on the 
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Fig. 1: Schematic of the experimental setup. LO pump beam focusing lens (f = 150 mm), PP- 
KTP periodically poled potassium titanium oxide phosphate crystal 1x2x11.5 mm^, BD beam 
dump, SLM spatial light modulator, LI and L2 two lens symmetric imaging arrangement (f = 
100 mm) to enhance the spectral resolution with a magnification of 1 :6, prism compressor com- 
posed of four N-SFll equilateral prisms, BF bandpass filter 4 mm BG18 glass, SPCM single 
photon counting module with a two lens (L3, L4) imaging system. The inset shows the mea- 
sured down-converted spectrum overlaid with a schematic illustration of the frequency-bins for 
a ququart. Each of the gray shaded areas represents a single bin whose amplitude and phase can 
be manipulated individually. 

SLM through 

d-l d-1 
j=Q j=0 

where lu^l and (py are controlled independently. Since in our experiment there is no spatial 
separation between idler and signal modes, we address each photon individually by assigning 
M^{uj) to the lower frequency part and M'*(cj) to the higher frequency part of the spectrum. 
Coincidences of the entangled photon pairs are detected within a time window of about 100 
fs through up-conversion in a second PPKTP crystal ll23l . The detection crystal is tempera- 
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ture stabilized at 35 °C to maximize the up-conversion rate. The recombined 532 nm photons 
are then imaged onto the active area of a single photon counting module (SPCM, ID Quan- 
tique, idlOO-50-uln). Since the entangled photons are detected in coincidences through up- 
conversion, we define the measured state \ip)s = dujT{uj)al{u)dl{—uj)\0)i\0)s where T{uj) 
describes the joint spectral amplitude of SPDC filtered by the phase matching properties of the 
detection crystal. The measured signal after shaping and the up-conversion stage then reads 



S 



dujT{uj)M'{uj)M'{-io) 



for a continuous wave pump and is equivalent to 

2 



s = \{x\i^) 



Cl 



1=0 



(6) 



for the direct product state \x) = (Ej^i^lj)^) (Ej'=o ^i'l^')^) l^)^'^^ ®- 
combination of the SLM together with an up-conversion coincidence detection therefore real- 
izes a projective measurement. Different quantum protocols can thus be implemented by the 
corresponding choice of 



4 Quantum state tomography of maximally entangled qutiits 

At first, maximally entangled states are generated by Procrustean filtering [|26l where we equate 
the Cj coefficients in Eq. (HI) by adjusting the amplitudes \uY\ of Eq. ([5]). The Procrustean 
filtering applies in two steps. First, we place the frequency-bins of Eq. (|3]) such that the states 
\j)i\j')s are well separated on the frequency axis and that no coincidence events are detected 
when j 7^ j'. By doing so, we implement the orthogonality of the frequency-bins. Second, 
considering the down-converted spectrum (Fig. [T}, it is obvious that the shape of the spectrum 
in combination with the bin width Acoj and its central frequency coj defines the amount of 
coincidences for the states with j = j'. It is then the \j)i\j)s state with the lowest coincidence 
rate which defines in Eq. (|5]) the amplitude scaling \uY\ of the others. We measure net count 
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rates of 43 Hz for a qubit, 13 Hz for a qutrit, and 6 Hz for a ququart. In an iterative procedure 
of scaling the amplitudes followed by measuring the \ j)i\j)s states, we balance the coincidence 
rates until they become equal within their errors assuming Poisson statistic. This filtering is a 
trade-off between maximal count rates and purity of the maximally entangled state and leads to 
the desired equation of the cj in Eq. (H)). 

Quantum-state tomography then allows us to retrieve the density matrix pd of these states 
by performing projective measurements [|29l |30 I . As a tomographically complete set of basis 
vectors we choose, inspired by [13J, single bin states \xji) = \ji)i\j'i)s and superposition of 
two bin states \Xa,n,j2) = Kb'i)^ + e"'\j2)i){\j[)s + e'^'ljs)^) with jf , = 0, d - 1 and 
Ji^ < with the relative phase a^'^ = 0, |. To obtain a positive semidefinite pd we employ a 
maximum likelihood estimation BTI . first applied to a qubit in [|29l and then extended to qudits 
up to (i = 8 in [13] . The reconstructed density matrices up to dimension d = A are shown in 
Figure |2l As a measure of how close the reconstructed state pd is to a maximally entangled 
state pme, we computed the fidelity Fd = Tr{^/ y/p^Pdy/p^) Ell- When pd is a maximally 
entangled state, the fidelity is equal to one. The computed fidelities of the reconstructed states 
are F2 = 0.928 ± 0.010, F3 = 0.855 ± 0.010, and F4 = 0.781 ± 0.018. We estimated the 
statistical 2cr-error with a Monte-Carlo method by randomly adding normally distributed count 
rate errors to each measurement outcome and recomputing the fidelity. With increasing dimen- 
sion we have to implement more frequency-bins within the same spectral range. Because of the 
finite spectral resolution of the setup, the increasing overlap between adjacent bins leads to a 
decrease of the fidelity. 
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Fig. 2: From left to right: Reconstructed density matrices of a qubit, qutrit, and a ququart, 
based on background-subtracted coincidence counts. Shown are the real and the imaginary 
parts. The small, residual imaginary values are due to remaining dispersion between the 
frequency-bins. 

5 Bell inequalities for non-maximally entangled quciits 

More generally, we obtain non-maximally entangled states by varying the amplitudes in Eq. |4l 
These states are then studied with regard to their non-local properties. The non-local properties 
of quantum correlations are usually measured by a Bell parameter / whose value is restricted 
under the local variables assumption. In [3], Collins et al. (CGLMP) generalized Bell inequali- 
ties to arbitrary c/-dimensional bipartite quantum systems by defining a dimensional dependent 
Bell parameter 1^. If correlations between two space-like separated systems can be explained 
through local realism, then Id ^ holds for all > 2. Despite a left-open locality loop- 
hole in our detection method, the violation of the precedent inequality indicates the existence 
of non-classical correlations due to entanglement. A counterintuitive property of the CGLMP 
inequality is that for dimensions d > 3 the inequality is more strongly violated by certain non- 
maximally entangled states than by maximally entangled states ||4l|5]|. In order to compare the 
sensitivity of the Bell parameters hil) and 13(7) to an entanglement parameter 7 G [0, 1], we 
consider the following bipartite qubit and qutrit states 

IV'(7))^'^ = -j=L={\0)a\0)b + i\1)a\1)b), (7) 
VI + 7^ 
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|V^(7))(=') = ^=L=(|0)a|0)b + 7|l)A|l)i. + |2)a|2)b), (8) 
v2 + 7 

where we associated the idler and signal photon with subsystem A and B. For each of the 
two subsystems individual measurement settings a,b E {1,2} and measurement outcomes 
Aa, Bh = 0, ... (i — 1 are assigned. According to CGLMP [3], the generalized Bell parame- 
ter in (i-dimensions is then defined as 

W2]-i / 2k \ 
h = 1-:^ {+[P{A^ = B^ + k) + P{Br=A2 + k + l) 



k=0 



d- 1 



and explicitly reads 



+P(A2 = B2 + k) + P{B2 = Ai + k)] 
-[P{Ai = Bi-k-l) + P{Bi =A2-k) 
+P{A2 = B2-k-l) + P{B2 = Ai-k-1)]} (9) 

/2 = /3 = +P{A, = Bi)+PiBi=A2 + l)+P{A2=B2) 

+P{B2 = Ai) - P{A, = B,-1)- P{B, = A2) 

-P{A2 = B2-1)- P{B2 = - 1) (10) 

for (i = 2, 3. For all d > 2, the inequality Id < 2 holds for local realistic theories. The 
probability that the outcome Aa differs from Bb by k modulo d is thus given by 

P{A, = Bb + k) = Y, P{Aa = J, Bb = (j + k) modd). (11) 

i=o 

Analogue, we have 

d-i 

P{Bb = Aa + k) = J2 Pi^a = (j + k) modd, Bb = j). (12) 

j=0 

The Bell parameter is described by a combination of projective measurements on the states 

\Xm,n) = |m);\|n)^ definedby 

1 "^^^ / 271 \ 

l"')^ = 7^ 5Z ('T-^^"" ^ ""7 l-^'^^' ^^^^ 
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withm,n = 0, . . . ,d—l and a specific choice of detection settings a i = 0,0:2 = 1/2, /3i = 1/4, 
and /32 = —1/4 for maximally entangled states f3]. Allowing for the states in Eq. (|7]) and 
Eq. dH]), the individual joint-probabilities become a function of 7. In accordance with Eq. Q, a 
measured coincidence signal is given by 

P{Aa = m,B, = n,7) cx | (x„^,n|^(7))^'f (15) 

with the projective states \xm,n) = |"^)^4|'^)b of Eqs. (I13I14I) . For a specific detection set- 
ting a, b each coincidence measurement has to be normalized according to P{Aa = m, Bh = 
n, 7) /Ma,b,'y with Ma,b,j = Yjm,n=o ^i^a = ^Ti, Bb = u, 7). Starting from a maximally entan- 
gled state through Procrustean filtering, the reduced entanglement, i.e. 7 < 1, is obtained by 
decreasing the transmission amplitudes of the bins associated with 1 1)^| 1)b using the SLM. 

It has been shown [5] for qutrits that the Bell inequality /s < 2 is maximally violated 
for -fmax = {VTl - V3)/2 ^ 0.792. We measured the Bell parameters Jf^ and /|^'^ for 
qubits and qutrits as a function of 7 (Fig. [3]). The experiment reveals a higher sensitivity to 
7 of the Bell parameter for qubits compared to qutrits, which is in accordance with theoret- 
ical predictions. The theoretical curves are scaled to the experimental data using the sym- 
metric noise model pXil) = -^d|V'(7))*-'^'' ^'^^(V^(7)l + (1 ~ Ad)1^2/(i^ in which deviations 
from a pure state due to white noise are quantified by a mixing parameter and 1^2 de- 
notes the (i^-dimensional identity operator. The value of the Bell parameter for p^"(7) then 
scales as Id{pT{l)) — ^dldil)- We experimentally determine the mixing parameters to be 
X'^^P = 0.920 ± 0.013 and X^^^ = 0.807 ± 0.008 where the 2(T-errors are calculated assuming 
Poisson statistics. Note, that the specific detection settings are not optimal for d = 2. In Figure 
[3l we therefore additionally depict values of the Bell parameter for optimal settings given by 
Horodecki's theorem [|33l . Similar to the measured 12(7) and in contrast to 13(7), the Horodecki 
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Fig. 3: l2^^ (blue dots) and J|^^ (red diamonds) show the experimental values of the Bell pa- 
rameter in dependence of the entanglement parameter 7. The 2cr-errors are calculated assuming 
Poisson statistics on background-subtracted coincidence counts. Straight lines show the the- 
oretically predicted Bell parameters 72(7) (blue), 12(7) using the Horodecki theorem (dashed 
blue), and 13(7) (red). The curves are scaled with their corresponding mixing parameter. The 
dashed black line indicates the local variable limit. 

curve is monotonically decreasing for 7 < 1. In the d = 3 case the optimal settings were only 
determined for 7 = 1 and 7max Recent numerical studies in [|34l report an even stronger 
violation of local realism by two entangled qutrits provided that a more general measurement 
bases is applied than the settings used in the CGLMP inequalities. 
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6 Conclusion 



By exploiting the flexibility of a SLM, we have been able to reconstruct maximally entangled 
qufiit states up to d = 4 through quantum state tomography. We moreover demonstrated the 
existence of frequency entanglement by measuring a CGLMP Bell parameter above the local 
variable limit for maximally and certain non-maximally entangled qubit and qutrit states. In 
our actual experiment, the available dimensions to encode qurfits in the frequency domain are 
currently constrained by the finite resolution of the optical setup and the pixel size of the SLM. 
The spectral resolution can be improved from 9 nm to 0.2 nm by replacing prisms with grat- 
ings such that the accessible dimension becomes then only limited by the number of pixel of 
the shaping device. This modification would result in a seven-fold increase in dimension. The 
here discussed method allows for an easy access to the spectral components of the photons and, 
further, a broad class of transfer functions can be implemented on the SLM. Therefore, other 
qufiit encoding schemes like time-bins and realizations based on Schmidt modes can be carried 
out. In addition, with shaped frequency entangled qudits, it is possible to investigate recently 
proposed theoretical results: Hilbert space dimensions can be probed using new dimension wit- 
nesses [|35l and extended Bell-tests with generalized entangled qutrit states are within reach 

m. 

Ultimately, the dimension of the Hilbert space is limited by the bandwidth of the pump laser 
and can reach values as large as a few millions. This provides a vast potential of encoding high- 
dimensional qufiits for quantum information and communication technologies. 
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